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We compute the spectral correlation functions for the transition from a harmonic oscillator 
towards the Gaussian Unitary Ensemble (GUE). We use a variant of the supersymmetry method 
to obtain analytical results in a fast and elegant way. In contrast to certain related transitions, 
the fe-point correlation function possesses the k x k determinant structure of the GUE limit for the 
entire transition. The results are used to consider also the spectral correlations of a superposition of 

] M transition spectra. Our results are non-perturbative and are valid for all values of the transition 

■ parameter. 
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I. INTRODUCTION 



Random Matrix Theory jl| is a powerful tool for the modeling of spectral fluctuation properties. Due to the general 
symmetry constraints, a time-reversal invariant system with conserved or broken rotation invariance is modeled by 
the Gaussian Orthogonal (GOE) or Symplectic Ensemble (GSE), respectively, while the Gaussian Unitary Ensemble 
' (GUE) models the fluctuation properties of a system under broken time-reversal invariance. These ensembles are 
known to describe the generic fluctuation properties of chaotic quantum systems very accurately [||-|| . While numer- 
t — ■ ical simulations of the ensuing matrix models usually pose no serious difficulties, the analytical calculations of the 
observables, i.e. the correlation functions, is generally a non-trivial task. In the case of the pure ensembles, Mehta 
and Dyson [jjj solved the problem about thirty years ago by introducing the orthogonal polynomial method. 

However, since a generic physical system has, classically, a mixed phase space, the spectral fluctuations of the 
corresponding quantum system will be in between the pure cases. The transition from preserved to broken time 
£>V reversal invariance was worked out by Mehta and Pandey ffl in 1983. The transitions of the spectral fluctuations in 
the case of gradually broken symmetries, i.e. quantum numbers was computed in Refs. [|| and 0. Here we will focus 
on a system that undergoes a transition from regular to chaotic fluctuations. To model such a system, we write the 
N x N random matrix representing the total Hamiltonian as a sum of a regular and a chaotic part, 

O ■ H(a) = H {0) + aH {1) (1) 



where a is the dimensionless transition parameter. The matrices H^- 1 ' are drawn from a Gaussian Ensemble with the 
probability density function P^\H^). Here, we are mainly interested in the transition from a regular, equispaced 



. N . 

spectrum to a chaotic one. This is a very important physical situation since many systems, particularly in nuclear and 
molecular physics can be described as a chaos producing part coupled to a harmonic oscillator. Often, it is necessary 
to give the regular part a block structure which reflects the presence of symmetries. This was already realized 
by Porter and Rosenzweig H who investigated, experimentally and numerically, atomic spectra that contain various 
angular momentum and spin quantum numbers subject to different coupling schemes. 

However, at the moment, we make no assumptions yet for the probability distribution pj^\H^) of the matrices 
H(°\ The decomposition (|l|) can be justified for potential and billiard systems. Detailed numerical simulations for 
the transition of the fluctuations can be found in Ref. ||. However, despite several attempts, full-fledged analytical 
discussions could only recently be performed for the case of broken time reversal invariance in which the matrices 
are drawn from the GUE. For a history of the studies devoted to these problems see Ref. [||. Presently, there 
are the following techniques which make such calculations possible. Pandey Q| presented a certain construction of the 
solution of Dyson's Brownian Motion Model H . A related approach was more recently put forward by Forrester jl0| . 
A very direct and compact technique for the GUE was constructed in Refs. [jll 12|. It relies on a variant of the 



supersymmetry method ]l3| , pT| which was introduced in Ref. The enormous simplifications are due to the fact 
that supersymmetry can, loosely speaking, be viewed as the "irreducible representation" of Random Matrix Theory 
which becomes apparent in a new class of diffusion equations Recently, Brezin and Hikami |l(| presented a 

third approach to derive similar integral representations. 



1 



In this paper, we will apply the methods of Refs. PIP^I to the transition starting from a harmonic oscillator. 
Pandey Q gave already a formula for the two-level correlation function on the unfolded scale. Forrester |l(J extended 
the result for the A;-level correlation function. Here, we have three goals. First, we will show that this result can be 
obtained very fast in a direct application of the general results of Refs. (ll] 121. Second, we will present plots and 
a detailed discussion of the two-level correlation function. Third, we will go beyond the known results and study a 
block structure of by considering the superposition of M transition spectra. After briefly sketching the method 
m Sec. g, we work out the crossover transition from one harmonic oscillator to the GUE in Sec. UI1 In Sec. nM 
study the superposition of M spectra. We summarize our results in Sec. |v|. 



we 



II. TRANSITION ENSEMBLES 



Before turning to the harmonic oscillator, we briefly summarize the general results. A detailed discussion can be 
found in Ref. Jl|], see also Ref. 0. 

As functions of the transition parameter a, we wish to study the fc-level correlation functions 



R k ( Xl ,...,x k ,a) = -1 fd[H^]P { N °\H^)x 



k 



d[HW]pg\HW) nimtr 1 (2) 
p= i x p \ a ) 



depending on k energies x p , p = 1, . . . , k. where the energies are given imaginary increments such that x p — x p ±ie. It 

is convenient to work with correlation functions R k {x±, . . . ,x k , ex) which involve the full Green functions, including real 
and imaginary parts. By studying different combinations of the signs of the imaginary parts of the Green functions, 
we can construct the physically interesting functions see Ref. The correlation functions R k {x\, . . . , x k , a) 

that may be written as the derivatives 



1 d k 

R k {x u ...,x k ,a) = TTTTk k Z k (x + J,a) 

\ ' i\. P =i aj p 



(3) 

j=o 



of the normalized generating functions Z k (x + J, a) . The energies and the source variables are ordered in the diagonal 
matrices x — diag^i, x\, ■ ■ ■ ,x k , Xk) and J = diag(— Ji, +J\, ■ . ■ , — Jfc, respectively. The desired functions 

Rk(xi, . . . , Xk, ol) can be derived from the generating function ^Z kl where the symbol 3 stands for the proper linear 
combinatio n K2| . The average over the GUE is done by means of the standard techniques of the supersymmetry 
method ]l3| , |14| , yielding 

Z k (x + J,a) = J d[HW]PP(H<®) J d[a] exp (--LtrgaA x 

detg- 1 ((x ± + J -a)®l N -l 2 k<8>H( A . (4) 

For details of the derivation and notation, the reader is referred to Ref. [{Isf . In Eq. (||), a denotes a 2k x 2k Hermitean 
supermatrix, and Iat and l2 k are N x N and 2k x 2k unit matrices, respectively. 

To proceed, it is in our case advantageous to avoid the saddle point approximation of Refs. Jl3|,|l4]]. We shift the 
matrix x + J is shifted from the graded determinant to the graded probability density and the supermatrix a is 
diagonalized according to a = u~ 1 su, where s = diag(sn, isi2, . . . , Sfei, is k 2). The volume element can be rewritten 
as d[a] — Bfc(s)d[s]d(J,(u) with B k (s) — det[l/(s p i — is g 2)]p,g=i,...,fc the Jacobian, here referred to as Berezinian. The 
non-trivial integration over the unitary diagonalizing supergroup with its Haar measure djJ,(u) is the crucial step 
and can be performed with the supersymmetric extension [fijl of the Harish-Chandra Itzykson Zuber integral p8| . 
Collecting everything we arrive at 

Z k (x + J,a) = 1 — rj(x + J) + 1 [ G k {s - x - J,a)Z<£\s)B k {s)d[s] , 

±> k {x + J) J 

Zf (x + J) = J d[H ]pW (IZo)detg- 1 ((x± + J) ® 1 N - l 2k ® H ) , (5) 
where the kernel resulting from the group integration is Gaussian and given by 
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G k (s -r,a)= 1 exp ( — ^trg(s-r) 2 ) . (6) 

with r = x + J. The distribution 1 — rj(x + J) in Eq. (^) ensures the normalization ZjAx, a) = 1 at J = 0. This 
distribution is not important for any of the formulae to follow, see the discussion in Rcf. |12]. The generating function 
Z k (x + J, a) satisfies an exact diffusion equation in the curved space of the eigenvalues of Hcrmitean supermatrices. 
Here, t = a 2 /2 is the diffusion time and the generating function Zy'ix + J) serves as the initial condition. This 
diffusion is the supersymmetric analogue |Q of Dyson's Brownian Motion Q . 

The integration over s requires to take a new type of boundary contributions []l3| , [l4 19 into account which do not 
occur in ordinary analysis. However, in Refs. |l5|,^(J it was shown that we do not need to worry about them when 
calculating correlation functions of the type we are interested in here. Collecting everything, we obtain the fc-level 
correlation functions 

Rk{xi, . . . ,Xk,a) = 

[ G k (s - x, a) SZ^ (s)B k (s)d[s] (7) 



for non-zero a. The case a = is trivial by construction. 

As a last step it remains to unfold the correlation functions for large level number N by removing the dependence on 
the level density. We define new energies £ p — x p /D, p — 1, . . . , k in units of the mean level spacing D. The transition 
parameter has to be unfolded, too, A = a/D and was first introduced by Pandey pl| . The fc-level correlation functions 
on the unfolded scale . . . , A) = limjv— »oo D k Rk{x\, . . . ,x k , a) are then generic, i.e. translation invariant over 

the spectrum. It is useful to unfold the integration variables s in Eq. <Yn) by making the rescaling s — » s/D. We arrive 
at 

Xk(€i, ■ ■ ■ ,£fe> A) = 

f G k (s - £, A) 34 0) (s)B k (s)d[s] (8) 

for non-zero A where the unfolded generating function of the arbitrary correlations is given by 

Z f( S )= hm Z^(Ds). (9) 

Hence, we have expressed the unfolded fc-level correlation function for the transition from arbitrary to GUE fluctuations 
as a 2fc-fold integral. 



III. TRANSITION FROM AN EQUISPACED SPECTRUM 

All results derived so far are correct for arbitrary initial correlations R k °\xi, . . . ,x k ) or JQ ' . . . We now 
apply them to the case of a equispaced spectrum. The harmonic oscillator probability density function reads 

4 0) (# (0) ) = 

N 

J] - e„) J] ^ReffW^ImffW) (10) 

n— 1 n>m 

where 

( (n- l)D/2 + SD, n = 1,3,5, ...,N 

e« = <^ (11) 
{ -nD/2 + 5D, n = 2,4,6,. ..,N- 1 . 

D = 2^/2/N is the mean level spacing and the number of levels, N, is odd. The spectrum is shifted by 6 and we 
require \6\ < 1. We find from Eq. (^) for the generating function of the initial condition 

7 (Q) ( s A isp2 -S {N yj /2 (l-(is p J(n + 6)){l + (is P2 /(n-6)) 
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According to Eq. (|9|) we evaluate the generating function on the unfolded scale in the limit N — > oo and obtain 

4 o )(s) = sin^s - 5) 
p = l sm7r(s Pl - S) 

The signs are determined by the choice of the sign of the imaginary increment in the Green function. The correlation 
functions can be worked out by using the general result (||) . The initial condition takes the form 

% z k ( s ) = TT sin n ( is P 2 - 5 ) Im T- FT • (14) 

p = i smTT{s pl -5) 

We use the identity 

Im - = V (-l) k S(s-k-S) (15) 

sinTr (s~ - 6 ) ^ 

K ' fc=-oo 

and insert Eqs. (^) and ( |l4|) into Eq. (Q). Since the function Bk(s) is a determinant and since the generating function 
of the harmonic oscillator spectrum ([l!]) and the Gaussian kernel (^) are products of 2k factors, the fc -point spectral 
correlation function may be written as a determinant |l5| 

X k fa, ■ . ■ , A) = det [C(&, A)] p>9=1 fe (16) 



with 



A) 



+ 00 + OC 



7T 2 A 2 J 7 s p i - is q2 
x exp 



(^ 9 2-^) 2 -(.S p l-^) 2 )) 

sin7r(is g2 - <5) 2J (— — / — <5) 

+ oo 



^ ^ v - exp A2 

I — — OO 

+ OG 

f sinir(is q2 -5) ( (s q2 + i£, q ) 2 \ , 

x J ds * ^TWTs) exp { a 2 J • (17) 

— oo 

The remaining integration may be done as follows. First, we translate the path of integration about +id which does 
not change the value of the integral. Then, we write the denominator of the integrand as a Laplace transform and 
perform the integrations by means of p2| 

, sinh7rs / (s + i(£ — 5)) \ 

ds — exp -5 = 

s + il I A 2 I 

— OO * ' 

-*(-!)' exp {^P^-) I-erf (^|^sign( - fy) . (18) 
This yields (erf denotes the error function as defined in Ref. [p2[) 

r( c c ^ 1 ( - & f - s ? 



^A 2 ^ V A 2 



x 
i=- 



exp ( -Tzitq-Zp) 1 



xlmerfj ^ '- sign(Z) - i-A ) . (19) 
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To proceed, we use the definition of the error function and write it as the integral of a Gaussian. Taking its imaginary 
part and writing the just introduced Gaussian as a Fourier transform yields 

C(£p; £gj A) = 

1 oo 

1 f , fn 2 X 2 t 2 \ f , 

at exp — - — / as exp 



4tt 



x 

i 



]T cos (tt(^ - 6 - l)t) cos - J - Z)s) . (20) 



We use Poisson's sum formula and perform the remaining integrations. The final result displays C^p, £ 9 , A) as a 
Fourier series 

oo 

C(£p,£ 9 ,A) = Y, exp (-tt 2 A 2 Z 2 ) 

/— — OO 

/ , , s , sinh(7r 2 A 2 / + i7r(e„ - , s 

xRe [e M ^ P - *)) ^ - _M j , (21) 

which inserted in Eq. ( [l6| ) completes our calculation. Two interesting limits may be considered. In the limit of 
vanishing transition parameter one obtains 



hmC(^,AH Sin ; ( ^ y £ 5(^-5-Z) , (22) 



which generates the spectral correlations of the harmonic oscillator. The opposite limit of infinite transition parameter 
yields the GUE spectral correlations, i.e. 



_ sinTrfe, -jg) 



lim C(£p, A) = -^^—1^ . (23) 

Plots of the level density 



AM£,A)=C(£,£,A)= £ exp(-7r 2 A 2 / 2 ) cos 2^- <5) (24) 

^ — — oo 

are shown for different values of the transition parameter A in Fig. [I]. Note that very small values of the transition 
parameter A already cause a considerable broadening of the (5-peaks of the original harmonic oscillator spectrum. It 
has been observed in almost all crossover transitions that, on scales of a few mean level spacings, the GUE limit is 
reached when A is of the order of unity, see the review in Ref. § . Figs. ^ and || show that this is also true for the 
two-point function 

X 2 (a,6,A)=X 1 (ei,A)X 1 (6,A)-C(ei,6 ! A)C(6,6,A) . (25) 

The delta-peaks at integer values £i ^ £2 are broadened for a small value of the transition parameter A, and the 
GUE-correlations are approached fast. Note also that the two-point function (|2^) becomes translation invariant (i.e. 
it depends only on the distance £1 — £2) if averaging over 6 is done p0|. 



IV. SUPERPOSITION OF SPECTRA 



We now consider a superposition of independent, non-interacting spectra as the initial condition. This is equivalent 
to saying that there is a symmetry that allows one to write the initial condition in the form — di&g(H[ \ . . . , H^) 
where each of the N m x N m matrices Hm\m — 1, ...,M is drawn from an independent ensemble. We have 
Yl m —i N m = N, the total dimension, and Yl m =i ^rn = where D m is the mean level spacing in the m th ensem- 

ble. The quantities g m = D/D m are referred to as fractional level densities. If all N m are equal, we have_o m = 1/M. 
Symmetries which lead to such models are frequently found in nuclear, atomic and molecular physics |@M. 
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For given fluctuations of every sub-ensemble with index to, the fluctuations properties of the superposition can 
be worked out in a rather straightforward way. In the case of the nearest neighbor spacing distribution, the result 
is given in the article by Porter and Rosenzweig f|. More general discussion may be found in Ref. Q and in the 
review p|. 

Here, we address the case that all matrices Hm are drawn from harmonic oscillators. In discussing the influence 
of a chaotic admixture to this initial condition, two different scenarios are of interest: The chaotic admixture can 
either preserve the symmetry or break it. In the first scenario, we have to add a GUE matrix to every matrix ■ 
Thus, we still have a superposition of independent spectra and the total matrix H(a) still has the block structure, 
i.e. H(a) = diag(ffi(a), . . .,H M (a)). 

In the second scenario, we add one N x N GUE matrix to the block diagonal matrix H 1 * ' . Thus, for non-zero 
transition parameter, the total matrix H(a) has no block structure. This scenario is in the spirit of the Porter- 
Rosenzweig model ||. For M = 2 such a symmetry breaking was investigated in Refs. g| and 0. In these studies, 
how ever, and in contrast to the present one, the initial c onditi on was also chaotic. We discuss the first scenario in 



Sec. IV A, before we briefly turn to the second one in Sec. IV B 



A. Symmetry preserving case 



Each of the block matrices is a harmonic oscillator, coupled to a GUE. For simplicity, we assume that all blocks 
matrices shall have the same dimension N m and the same coupling to the GUE. Moreover, to lift the degeneracies, 
the m th oscillator spectrum is chosen with an energy shift 5 m — m/M (cf. Eq.(|ll|)). We are interested in the unfolded 
two-point correlation function 



M 



X 2 (&>£2,A) = ^ 9m X 2,m{9mix,gm^2,9m>) 



m—1 
M 



+ J]] gm9nXi t m(g m £i,g m X)Xi >n (g n {;2,gnX) 
n,m— 1 
M 

~ 9rn.Xl,m(gm£l, ffm ^)-^l,m (gm^2 ; 9mty ■ 



(26) 



where g m = X/M is the fractional density of each block matrix. To motivate this equation, we notice that 

X% (^1,^2, A) is the probability density of finding a level at £i and another one at £2- The levels can either be- 
long to the same symmetry sector (first line of Eq. (|2^)) or to different symmetry sectors (second and third line of 
Eq. (^6|)). Since the second line includes also correlations within one symmetry sector (these have been accounted 
for in the first line), the third line must be subtracted. We emphasize once more that the functions Xi tm (g m ^i, g m X) 
are not unity in the case of a harmonic oscillator. 

We now insert the the two-point function ( p5|) with the appropriate energy shifts into Eq. ( 26 ) and interchange the 
summation over m with the summation over k, I. For M > 1, this yields expressions like 



A I 



cos2 7 r-^-(£ 1 - to) cos 27r-^- (£> -m) = — (S, 



k,l 



5 k _i)cos27r— (£1 -£ 2 ) 



(27) 



which help to perform one more summation. The final result reads for M > 1 



sinh (n 2 X 2 k/M 2 + iirfa - &)/M) 



Tr 2 X 2 k/M 2 + i7r(£i — ^)/M 



o.xp j -2tt 2 A 2 -^ ) cos2^ — I'., - < ; ,i 



In the limit of vanishing A we find 



hm4 M) (a,6,A) = i-5(a-6) 



M 



(28) 



(29) 
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and in the limit of infinite transition parameter we recover the two-point correlation function for a superposition of 
M GUE's 



a->oo (a,6,A) = i--(^ <ix _^ m ) ■ (30) 

For large M and nonzero transition parameter A, the two-point correlation function approaches the value one. 
Thus, the spectrum becomes completely uncorrelated and exhibits the same correlations as the (randomly gener- 
ated) Poisson-spectrum. We also note, that the two-point correlation function ( psj ) depends only on the difference 
(£1 — £2) of its arguments. This translation invariance results from the assumption that the m th spectrum is shifted by 
m/M mean level spacings. Thus, any deviation from translation invariance gives a hint to clustering or degeneracy of 

levels belonging to different symmetry sectors. Figure [I] shows the a plot of M[l — X^^i, £2, A)] for various values 
of the coupling A. This function depends only on the ratios X/M and (£1 — t^/M. 



B. Symmetry breaking case 

We turn to the second scenario defined above. We do this briefly, because we mainly want to show the very different 
structure of the result. In general, the initial condition on the original scale takes the form 

M 

TCI— 1 

where the initial conditions zf^is) are still completely arbitrary. On the unfolded scale, we have 

M 

£\*) = U z ki(9ms). (32) 

m— 1 

These formulae have to be inserted into Eqs. (Q) and (g), respectively, to obtain the correlation functions. 
If all M initial ensembles are chosen as harmonic oscillators, we find from Eq. ( |l3| ) on the unfolded scale 

4 o) (*) = nn ^'M - (33) 



Since this is still a product, the determinant structure (|16|) of the fc-point correlation functions is still preserved. The 
functions 



C(£p> £q, A) 



+00 +00 
1 f f ds„i ds n 



7T 2 A 2 



s p i - is q2 



— 00 —00 



x exp ^ ((is q2 - - (s p i - £p) 2 )^) 

M M 



TT simr(g m is q 2 - 5 m ) Im TT = — - . (34) 

Sln7r (f™' s 9 l -°m') 



X 

m— 1 



are the entries of this determinant. For general g m and 8 m , this result is not easily amenable to simplifications. In 
some special cases, however, the product over m in Eq. (B3[) can be performed. In particular, for g m = 1/M and 
5 m = m/M we obtain 



4 o) w = n^% (35) 



. Sill 7TS 

P =i 



by using standard results for trigonometric functions. This is precisely Eq. (|13|) for 6 = 0. Thus, we recover the 
results of Sec. Ill because the case discussed here corresponds to the coupling of one oscillator, with the spacing on 
the original scale multiplied by 1/M, to one GUE-matrix. 
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V. SUMMARY 



We have computed all spectral correlation functions for the transition from a harmonic oscillator to the GUE. 
We have used these results to compute the two-point correlation function for a system which undergoes a crossover 
transition to a GUE starting from a superposition of independent oscillator spectra. We discussed chaotic admixtures 
to an initial block structure which either preserve or break this symmetry. These ideas and techniques can also be 
used for other block diagonal initial conditions with different individual statistics. 

Applying a variant of the supersymmetry technique, we gave an elegant and very straightforward derivation of these 
results. 
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FIG. 1. Level density A) for different values of the transition parameter A 
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X 2 (£i,&.A= 0.25) 




FIG. 2. Contour plot of X 2 (£i,£ 2 , A = 0.25) 
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X 2 (£i,&,A= 1.0) 




FIG. 3. Contour plot of X 2 (£i, 6, A = 1.0) 
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FIG. 4. Plot of M[l — X 2 (£1,^2, A)] for different values of the transition parameter A. 
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